Fractal dimension and threshold properties in a spatially correlated percolation model 



o 
o 

(N 
< 



C/3 



I 

O 

o 



(N 
> 

O 



00 
O 
On 
O 



X 



Hongting Yang/'^ Wen Zhang/ Noah Bray- Ah/' ^ and Stephan Haas^ 

'Department of Physics and Astronomy, University of Southern California, Los Angeles, CA 90089-0484 

'^Department of Physical Science and Technology, 
Wuhan University of Technology, Wuhan 430070, P.R. China 
^Department of Physics and Astronomy, University of Kentucky, Lexington, KY 40506 

(Dated: August 9, 2009) 

We consider the effects of spatial correlations in a two-dimensional site percolation model. By 
generalizing the Newman-Ziff Monte Carlo algorithm to include spatial correlations, percolation 
thresholds and fractal dimensions of percolation clusters are obtained. For a wide range of spatial 
correlations, the percolation threshold differs little from the uncorrelated result. In contrast, the 
fractal dimension differs sharply from the uncorrelated result for almost all types of correlation 
studied. We interpret these results in the framework of long-range correlated percolation. 
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I. INTRODUCTION 

Percolation refers to the formation of long-range con- 
nectivity across a lattice or network, and applies to 
various complex physical systems, including aerog elfll, 
polymers [2], superconductors [5, 0, IE Hi j and magnets [7, 
H]- Outside physics, percolation also arises in the context 
of ecology [p, complex networks [l^, [ll|, [H, [l^l and social 
science 

In these applications, sites or bonds are assumed to 
randomly occupy an underlying lattice or network. At 
a critical concentration, the percolation threshold p^, a 
connected cluster of sites spans the system. If the oc- 
cupation probabilities on different sites are independent, 
the problem maps onto the g-state Potts model in the 
limit q — » In two dimensions, the critical proper- 

ties of the g-state Potts model can be obtained from an 
equivalent Coulomb gas system for arbitrary In the 

limit q — > 1, one obtains the percolation length exponent 

= I and the fractal dimension = || of the percola- 
tion cluster appearing at threshold. A rigorous derivation 
of the same results follows from a map of percolation onto 
Schramm-Loewner evolution with k = 6[l9|. 

In physical applications, it is common for there to be 
correlations between the occupation probabilities on dif- 
ferent sites. This case can also be mapped to a g-state 
Potts modelf^O^I; however, the resulting model has corre- 
lated coupling-constant disorder. For short-range corre- 
lations, the disorder is irrelevant and does not change 
universal properties such as ly and Df. Interestingly, 
this is also true in the case where the correlations de- 
cay asymptotically as a:~° for sufficiently large distances 
X, provided that the exponent satisfies a > D where D 
is the number of spatial dimensions [20|. 

For long-range correlations, the critical properties may 
change. In particular, for power-law correlations satisfy- 
ing the extended Harris criterion av — 2 < 0, a long- 
range correlated renormalization-group fixed point has 
been predicted to occur with new percolation length ex- 



ponent l^lo 



and new fractal dimension [2C 



Al- 



though estimates have been made of the fractal dimen- 



sion near the upper critical dimension I? = 6, no general 
predictions or numerical work exist in two dimensions, 
the lower critical dimension. This work provides one of 
the first estimates of the fractal dimension in the presence 
of long-range correlations in two dimensions. Prior work 
has evaluated the fractal dimension of critical Ising clus- 
ters on the triangular lattice, a particular case of long- 
range correlated site percolation with a = ^ [131 • 

We focus on spatial correlations in a site percolation 
model on a square lattice. There are various ways of 
implemen ting correlations. For example, in an earlier 
study[2TI, I22I, a vacant site was assigned a probability 
R for being occupied. Here is a number representing 
enhancement, R is the probability for occupation for each 
of those vacant sites not having an occupied neighbor, 
and y is the number of its occupied neighbors. F = 1 
corresponds to the uncorrelated percolation situation; for 
-F > 1, the vacant sites with more occupied neighbors 
have larger occupation probabilities than the others. 

This model can be used to simulate correlated absorp- 
tion. An interesting result found for this model was that 
the percolation threshold pc is a non-monotonic function 
of the correlation strength F. In particular, correlations 
between sites do not always promote the occurrence of 
percolation, but in some cases instead shift the percola- 
tion threshold to larger occupation densities compared to 
the uncorrelated case. 

To investigate long-range correlations, a more efficient 
algorithm for growing the percolation clusters is needed. 
A few years ago, Newman and Ziff developed a fast Monte 
Carlo algorithm for site and bond percolation [2^. In this 
algorithm, a sample state with n -\- 1 occupied sites or 
bonds is realized by adding one extra randomly chosen 
site or bond to a sample state with n sites or bonds. 
In this way, a relatively small amount of computational 
effort is needed to determine the new clusters from the 
old. Here, we generalize the Newman-Ziff algorithm to 
the case of correlated percolation. 
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II. GROWTH PROCESS 

For simplicity, let us define the distance between near- 
est neighbor sites to be one unit and introduce an ad- 
justable integer length parameter d. In each configura- 
tion, starting from an empty square lattice, only the first 
occupied site is completely randomly chosen. The sec- 
ond site to be occupied is determined by first randomly 
choosing a vacant site, and then checking the projections 
of its distance to the first occupied site both along the 
X-axis and y-axis. 

If both projections are no larger than d, we occupy 
the vacant site. Otherwise, we randomly choose another 
vacant site and check its distance from the initial site, 
until finally both projections are smaller than or equal to 
d. In a similar way we occupy the other vacant sites. The 
correlated distribution is thus obtained by ensuring that 
there exists at least one occupied site with its distance 
to the next occupied site no larger than d, both in the 
horizontal and the vertical direction. A main difference of 
this correlated site percolation model to earlier ones[2ll 
[22j is that we do not permit the occupation of any vacant 
sites farther than d away from the previously occupied 
sites. 

In a square lattice of linear dimension L, d ^ L — 
1 corresponds exactly to the uncorrelated case, whereas 
d = 1 corresponds to the most compact population. For 
d = 1, a snapshot of a realization of ten occupied sites 
on a square lattice of L = 10 is shown in Fig. 1(a) which 
displays obvious spatial correlations. For comparison, a 
randomly populated realization for the same density is 



also shown in Fig. 1(b) 



Although the percolation process contains a length d, 
the correlations between site-occupation probabilities can 
in principle be long-ranged. We define site-occupation 
variables 9i at the sites {«} of a square lattice. The 6i 
take on the values 1 and corresponding to occupied and 
vacant sites, respectively. Averaging over cluster realiza- 
tions gives the site-occupation probability p = which 
is the same on all sites in the growth process considered 
here. To assess correlations, it is natural to compute the 
site-occupation correlation function: 



(1) 



From the homogeneity of the percolation process, it fol- 
lows that the strength of correlations between two sites 
i and j depends only on the displacement — ri — rj ■ 
Now, for a finite-size system with N sites, it is inter- 
esting to consider the case of site-occupation probability 
Pn — ct/N, where, the positive integer a <^ N does 
not depend on the system size N. Then, the correlation 
function ge(fij) ~p'n ~ —c?!^"^ tends to a non-zero 
constant for sufficiently large displacements ad. 
Physically, for such a small site-occupation probability, 
only a few sites a — NpN are occupied before the growth 
process stops. However, if sites i and j are sufficiently 
far apart such that ^ ad, then at least r^jd ^ a 
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(a) Correlated 
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(b) Random 

FIG. 1: (a) Realization often occupied sites in a 10 x 10 lattice 
for d = 1 (strongly correlated case); (b) the counterpart of a 
randomly populated lattice at the same density. The occupied 
sites are marked by "X" , the empty sites by "O" . 



steps of the growth process are required for both sites to 
be occupied. 

In any realization, then, it is impossible for both sites 
to be occupied, once the sites are sufficiently far apart. 
As a result, the product of site-occupation variables van- 
ishes, QiQj = 0, and the correlation function tends to a 



negative, non-zero constant geirij 3> ad) 



This 



shows that, at least for finite systems, correlations be- 
tween the site-occupation probabilities can extend over a 
long range, on the order of the system size. Turning to 
the limit N ^ oo, we expect that the correlation function 
ge [f'ij ) will decay to zero as oo for any p and 

d. In particular in the vicinity of percolation threshold, 
we anticipate that the correlations decay algebraically 
ge - l/r^ 



III. RESULTS 



The results for percolation threshold (Fig. 2(a)) were 
obtained by averaging over 3000 realizations for the 
largest L = 200 system sizes [13]. For a given realiza- 
tion, we find the occupation fraction Pc{L) when a span- 
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ning cluster first appears. Averaging over different real- 
izations, we arrive at a finite-size estimate that we an- 
ticipate has the following scaling behavior as L ^ oo: 
Pc{L) ~ Pc(oo) — bL~^/'^, where, v is the percolation 
length exponent [H, |24[. For uncorrelated percolation, 
one has the exact result v = ^ allowing a simple extrap- 
olation to infinite system size. For the correlated case, 
we extract the values of v and pc{oo) simultaneously. 

We begin by assuming v — 1 (See Fig. 2(a) ) and find 
a preliminary estimate of Pc{oo) « 0.793 ± 0.002, for the 
d — \ data shown in the figure, using a linear least- 
squares fit to the finite-size scaling form. To check for 
self-consistency, we take Pc{oo) in the equivalent finite- 
size scaling ansatz log(pc(oo) —pc{L)) = — Mog L -I- log &, 
and use linear least-squares fitting to extract v. The 
results for pc{oo) for various d shown in Fig 2(b) were 
obtained using v = \. 

The value of the percolation threshold for d ~ 1, 
Pc(oo) = 0.793 ± 0.002, is significantly larger than the 
uncorrelated value Pc(oo) = 0.5927[23|. As d increases, 
Pc{oo) decreases quickly, and reaches a minimum around 
d = 9, i.e. pc(oo) = 0.5910. 

It is natural that at small d a larger occupation fraction 
is required for percolation to occur. Percolation implies 
that distant sites are not only simultaneously occupied, 
but also joined in the same cluster. Reducing d, we find 
that distant sites have stronger anti-correlation of their 
occupation probabilities: they are less likely to be simul- 
taneously occupied. 

It is noteworthy that the percolation threshold pc{oo) 
in Fig. 2(b) gently increases with d beyond d = 9. For 
large d, an increasing number of the randomly added 
sites are not part of the spanning cluster, and therefore 
percolation requires a larger population density. Numer- 
ically we find that for d = 20, the percolation threshold 
is slightly smaller Pc{oo) = 0.5924 than the threshold at 
d = 40, Pc(oo) ~ 0.5927. The latter is indistinguishable 
from the uncorrelated result, at this level of precision. 
It would be interesting to consider other correlated per- 
colation processes for which this dip in pc{oo) is more 
pronounced. 

Remarkably, the extracted values of Pc(oo) depend 
very weakly on v, allowing for a precise estimate despite 
the large uncertainty in v: We find self-consistent fits 
using V = 1.0 ± 0.3 for all d. Further work is necessary 
to provide a reliable estimate of the correlated percola- 
tion length exponent v. Weintraub [l^l has suggested 
the result viong — 2/a, where, the exponent a governs 
the decay of correlations gg ~ 1/r" at large distance. It 
would be interesting to compare more precise estimates 
of V to this prediction, assuming that the decay of the 
correlations takes this form. 

To gain a better understanding of the effect of correla- 
tions near threshold, we turn to the fractal dimension Df^ 
governing the growth of the number of sites in the span- 
ning cluster s{oo) cx L^^ In Fig. 3(a) we plot the 
size of the spanning cluster on a log-log scale as a function 
of system size at fixed correlation d=l. The slope of the 
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FIG. 2: (a) The finite lattice percolation threshold Pc{L) 
varies with the inverse of the linear dimension l/L^''" for 
d — 1, with u — 1.0. Linear extrapolation of the values of 
Pc{L) at large L yields pc(oo) = 0.793 ± 0.002 at L = oo. (b) 
Dependence of the percolation threshold Pc{oo) on the corre- 
lation parameter d. With increasing d, pc decreases quickly at 
first, reaches a minimum at about d = 9, then grows slowly. 



line gives Df = 1.99 ± 0.04 for d = 1, indistinguishable 
from the Gaussian fixed point result Df = 2. Similarly, 
we obtain the fractal dimensions corresponding to other 
correlation strengths. The dependence of the fractal di- 
mension on the correlation parameter d is illustrated in 
Fig. 3(b) This plot shows a dependence similar to that 
of correlated percolation thresholds in Fig. 2(b). Increas- 
ing d from 1 to 9, the fractal dimension of the spanning 
cluster decreases, until a shallow minimum around d = 9, 
after which the fractal dimension increases towards the 
uncorrelated value. We note however a quantitative dif- 
ference: the fractal dimension Df — 1.80 ± 0.02 for the 
largest d = 40 data shown is still far less than the exactly 
known value for uncorrelated percolation 91/48 w 1.896. 
This is to be contrasted with the numerically obtained 
percolation threshold at d = 40, which coincides with 
that of uncorrelated percolation, within the precision of 
our numerical calculations. 
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FIG. 3: (a) The logarithm of the spanning cluster size at 
the correlated percolation threshold for d = 1, plotted as a 
function of the logarithm of the linear dimension. A linear 
fit of this plot gives the slope 1.99 ± 0.04. (b) The fractal 
dimension D, as a function of the correlation parameter d. 
The fitting errors are also given. 



IV. DISCUSSION 

We interpret the results for the fractal dimensions as 
an indication for a crossover, as function of d, between 
different renormalization-group fixed points. The frac- 
tal dimension Df = D/2 + 1 — rj/2, can be expressed 
in terms of the anomalous dimension rj governing the 
probability that two sites belong to the same cluster. 
Close to the upper critical dimension D = 6, three fixed 
points dominate the phase diagram for long-range corre- 
lated percolation [2^. For the Gaussian fixed point, the 
anomalous exponent rigauss = vanishes. For the pure 
fixed point (at which correlations between couplings in 
the Potts model are irrelevant) and the long-range fixed 



point, r] is non-zero. In contrast to the pure and Gaus- 
sian fixed points, the long-range value depends on a: 
Viong = (I? — 2 — a)/ 11 for a 4. Above the upper crit- 
ical dimension D > 6, the Gaussian fixed point governs 



for 



gauss 



2 > 0. Below the upper critical dimension 



D < 6, the pure fixed point governs once avp 



2 > 0. 



The long-range fixed point governs when neither the pure 
nor the Gaussian do. Assuming these results extend qual- 
itatively to the lower critical dimension D ^ 2, one may 
interpret our numerical results as follows. 

For d ^ 10, we observe Df = 2, corresponding to a 
vanishing anomalous dimension. This is naturally inter- 
preted as a Gaussian fixed point. In the limit d = 00, our 
growth process is identical to uncorrelated percolation 
and is governed by the pure fixed point with Df = 
In between these two extremes, for large but finite d, our 
numerical results are inconsistent with both the pure and 
the Gaussian fixed points, and are likely to be governed 
by a long-ranged fixed point with fractal dimension that 
depends on d through the exponent a. Within this pic- 
ture, the non-monotonic dependence of Df on d results 
from the presence this new, long-ranged fixed point at 
intermediate d. 

In conclusion, we have generalized the Monte Carlo 
method by Newman and Ziff using a simple correlated 
percolation model. We have then calculated percola- 
tion thresholds, spanning cluster sizes and fractal dimen- 
sions for different correlation ranges. The observed non- 
monotonic dependence of the percolation threshold on 
the correlation range can be understood in terms of two 
competing mechanisms: compactification of the spanning 
cluster as d ^ 1 and increasing abundance of particles 
that are not part of the spanning cluster as d — > cx3. The 
dependence of the fractal dimension Df on is less triv- 
ial. We have presented a scenario involving a cross-over, 
as a function of d, between three renormalization group 
fixed points. It would be interesting to apply this analy- 
sis to physical systems with tunable spatial correlations 
in the site-occupation probabilities, using it to control 
and optimize their critical properties. 
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